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THE EXPECTED DURATION OF RANDOM 
SEQUENTIAL ADSORPTION 
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When gas molecules bind to a surface they may do so in such a 
way that the adsorption of one molecule inhibits the arrival of others. 
We consider random sequential adsorption in which the empty sites 
of a graph are irreversibly occupied in random order by a variety of 
types of "particles." In a finite region the process terminates when 
no more particles can arrive. A universal asymptotic formula for the 
mean duration is given. 

1. Introduction. In this paper we consider a process known as random 
sequential adsorption (RSA) in wliich "particles" are deposited on a graph 
successively, subject to a blocking effect, whereby the presence of a particle 
at a particular set of sites may prevent the arrival of other particles at the 
same or nearby sites. The principal motivation comes from the adsorption 
of gas molecules onto a surface. Initially the graph is empty and, at rate 1, 
independently of all other particles and sites, each type of particle attempts 
to occupy each suitable set of sites. No particles leave the graph, so that 
once the arrival of a particle is blocked it is blocked forever. Thus, the 
process terminates on any finite set when no further particles can arrive. 
The most common types of RSA which have been considered are those 
in which single particles (monomers) arrive at single sites, but inhibit the 
arrival of neighboring particles, or when pairs of particles (dimers) can only 
arrive where pairs of positions are free. There are fairly general treatments 
of these processes in Penrose and Sudbury [2], a review article by Evans [1] 
and more recently the series of surveys in Volume 165 of the journal Colloids 
and Surfaces (2000). 

We wish to consider RSA processes on Z*^ which are sufficiently flexible 
to allow a variety of possible changes at any site. Thus, we do not restrict 
ourselves to monomers and dimers or to any particular ways in which the 
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arrival of one type may prevent the arrival of another type. To emphasize 
that there is no restriction on the shape of a particle, the term "config" is 
used in this paper. A configuration Tj is defined by a finite set of sites which 
must include the origin as "anchor." Any translation of this set is called a 
"config" of type Ti. A config arrives at a site x when its anchor attempts to 
arrive at x. If no neighbors of the config have arrived beforehand, then the 
arrival will be successful. We prove the following: 



Theorem 1. Suppose we have an RSA process in which k different types 
of config attempt to arrive at each site, each independently at rate 1. Con- 
sider a subset A of Z'^ consisting of n sites. Call the expected time until the 
process terminates on the n sites i5^term(^)- Call the set of subsets of Z'^ of 
n sites Sn- Define p to be the limiting probability as t ^ oo that an arrival 
at a site will be successful. If p> 0, then, as oo. 



max 



l + ^ + ^ + - + + Hp) - EtcrraiA) 



0. 



Exact values of p are not generally known, but in 1-d the dimer process 
has p = and the annihilation process p = e~^ (see, e.g., Sudbury [3]). In 
the case of dimers arriving on Z, it would be the limiting probability that, 
given that a dimer had not arrived at the pair of sites 0, 1, dimers already 
occupied pairs —2, —1 and 2, 3. Simulations of the dimer process on a set 
of 400 sites gave a mean that differed from the asymptotic estimate of 4.567 
by approximately 0.02. With this number of sites it can be seen that the 
ln(p) = — 2 term plays an important role. It will be clear that the variance 
of the expected time is bounded hy J2^/^'^ = 7r^/6. 



2. The expected time for the process. If there was no interference be- 
tween configs, then the expected time until all had arrived would be 1 + 
1/2 + • — h l/{kn). If the probability of the successful arrival of a config was 
p independent of all other configs, then the expected time until the last ar- 
rival would be J2r=o^ Pi^ ~ pY {J2i=r+i'''~^} which, asymptotically, tends to 
1 + 1/2 + • — h l/{kn) + ln(p) , as will be shown at the end of this section. The 
chief idea of the proof given here is that a subset of sites of 0(n) do behave 
as if they are almost independent when n is large. This is because, although 
the number of neighbors at a distance / in Z'^ grows as l'^, the infiuence of 
one site on another falls off as 1/V.. Further, the "bulk of the action" in the 
series 1 + 1/2 + • • • + l/[kn) occurs at the beginning, that is, for the later 
stages of the RSA process. 

It should be noted that the bounds given in the following lemmas do not 
depend on the positions of a set of configs, only on their number. 

Assume we have a finite set of types of finite configurations of sites 
{Ti, . . . , Tfc}. [E.g., in 2-d we might have the two configurations Ti = 
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{(0,0), (0,1), (1,1)}, r2 = {(0,0),(2,0)}.] The origin is chosen to be the an- 
chor. Any translation of a member of the set of types is a "config" [e.g., the 
set of points (m — 1, n — 1), (m — 1, n), (m, n) would be a config of type Ti]. 
The anchor translates with the config. For each config C there is a single 
arrival time t{C) which is a mean 1 exponential random variable, which is 
independent of the arrival times of all other configs. We may thus construct 
the process via a Harris diagram in the usual way. 

Two configs Ci, C2 are said to be neighbors if the successful arrival of one 
of them must block the successful arrival of the other. If a config has not 
arrived and one of its neighbors has made a successful arrival, then it enters 
a blocking mode. In a finite region, once all arrivals have taken place or all 
unarrived configs are in a blocking mode, the process terminates. 

Assume that initially all the sites of Z"^ are unoccupied so that the process 
is translationally invariant. Suppose the process on Z'^ has been running for 
time t. Define Pi{t) to be the probability that a particular config of type Tj 
is not blocked at time t given that it has not arrived. Since once any config 
is in blocking mode, it remains so, this probability is monotonic decreasing 
and limt^ooPiit) =Pi- 

In a realization, a config Cq is said to "affect" another config C, if there 
is a sequence of configs Co,Ci,C2, . . . ,Ci = C such that t{Ci) < t{Cj), i < j, 
and Ci is a neighbor of Cj+i. (Note that Cq "affecting" C does not imply 
that the arrival of Cq either ensures the blocking of C or its successful 
arrival. However, Cq "affecting" C is a necessary condition for the arrival 
or nonarrival of Cq to affect C in the sense that a realization with arrivals 
identical to the first, except for the nonarrival of Cq before C, will differ from 
the first in the successful arrival or nonarrival of C.) All distances are Li. 
Let Df, be the largest distance from the anchor of one config to the anchor 
of a neighbor. An m-neighborhood of a config is the set of configs whose 
anchors are distance < m from the config's own anchor. We put N = nk. 

Lemma 2. Define Am{C) to be the event that a config C is affected by 
at least one config outside its m-neighborhood. Then P{Ara{C)) < A^~™f for 
any positive integer niQ when m = [DbN^^"^^] and N is large enough. 

Proof. The probability that in a config path of length / the events are 
ordered in time is =l/l\. Let Nb be the maximum number of neighbors of 
any config, then the number of paths of length / is < (iVj)'. To be distance 
m away, the path must be at least length m/Dh. Thus, the probability a 
config is affected by another at distance m or greater is 



P{A^{C))< ^<2 




2d-| Afi/2<*/2d 



{[m/D,])l 



< 



N 



k=lm/Dt] 
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for any positive integer rriQ when m = [DbN^^'^'^] and N is large enough. □ 

If there are no config paths from outside the m-neighborhood, it is clear 
that events outside the neighborhood cannot influence the config, thus, we 
have the following: 

Lemma 3. Let C be the event that the config C arrives successfully. Let 
E he an event defined by sites outside the m-neighborhood of C. It follows 
that P{C\E,A<;^iC)} = P{C\A<;^iC)}. 

The number of configs in the m-neighborhood is called N^- 

Lemma 4. In d- dimensions, Nm < My/N where m = [DbN^/'^'^] and M 
is a fixed constant. 

Proof. The set of sites with Li distance from the origin < m is con- 
tained in the d-dimensional cube with each co-ordinate lying in [— m,m]. 
The number of sites is thus less than (2m -|- l)*^ and the number of anchors 
at each site is k. Putting m = [Di,N^/'^'^], the results follows. □ 

Let Qn he a, set of configs. In any realization of the RSA process a 
subset of r configs will be the last of the N to arrive. Define the event 

E{r,Qi\i,m) = {Each of the last r configs to arrive is 

the last of its m-neighborhood}. 

Note that arrival here does not imply successful arrival. 

Lemma 5. If N is large enough so that rNm < N, then 

r'^N. 

P{E{r, QM,m)}>l- > 1 - M-^. 

Proof. Call the last r configs to arrive Ci,C2, . . . ,Cr. There are at 
least configs available to be Cr, and one of them must be the last of 
at most + rNm configs. Cr arrives after Ci, C2, . . . , Cr-i, and since each 
Ci is also the last of its m-neighborhood, C,. must arrive after the union of 
these m-neighborhoods. Thus, the probability of this occurring is > N/{N + 
rNjn)- Conditional on this event, we now calculate the lower bound of the 
probability for Cr-i to be the last of its m-neighborhood. There are at least 
N — Nm of the original available to be Cr-i and one of them must be 
the last of at most N + (r — l)Nm configs. Thus, the conditional probability 
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that the second last is the last of its m-neighborhood is > (A^ — Nm)/{N + 
(r — l)Nm)- Continuing the argument, we obtain 

P{E{r, QN,m)} > 



{N + riV„)(iV + (r - l)N^) ■ ■ ■ (N + N^) ' 
Now, 

N-kNm rNm „ / rN„ 



> 1 - ^ ^ P{E{r,Q^,m)} > (l " 



N + {r- k)Nrr 
from which the result follows using Lemma 4. □ 

Define p = Ylj=iPj/k- 

Lemma 6. Given none of the last r + 1 configs to arrive is in the m- 
neighborhood of any other, then the probability all except the first of them 
are blocked = B^'^^p{l — pY , where ""^ < B < e^ . 

Proof. First, note that the assumptions imply that each of the last 
r + 1 configs to arrive is the last of its m-neighborhood. Call the last r + 1 
configs Ci, C2, . . . , Cr+i and define Ei, z = 1, . . . , r + 1, to be the event that 
Ci is blocked. We now show that the dependence of each event on the others 
is small, because being affected by events outside the m-neighborhood is 
small: 



p[E,\f]Ej/E, 



Pi E,. 



An(C,), fl E.j/Ei P A,n{C{) I fl Ej/E, 



+ P\Ei 



ai^{c,),{^e,/eAp[ai^{cA{^e,/e, 



p[e, An(C,),n^j/^*V(^m(Ci)) + P(i?i|^^(C*))^'(A 



m)' 



using Lemma 3 and the fact that the event AmiCi) is determined by events 
in the m-neighborhood of Ci. The probability of Ci being blocked uncon- 
ditioned on the status of the other r configs differs only in the first term 
above, which is less than P{Am)- Thus, putting the conditional probability 
of being blocked as + (1 — Pi), where may be negative, and averaging 
over the {pi}, we obtain e+ (1 — p), |e| < P{Am)- This equals B{l—p) where, 
using Lemma 2, we have e~^ ""^ < S < . □ 
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The values of B will depend on time and on the positions of the configu- 
rations. In the interests of simplicity of presentation, the values of B are not 
always the same in what follows, but they all satisfy the above inequality. 

Lemma 7. The probability the last r arrivals are blocked is < (3'^ for fixed 
13 <1 when r < kn/2Nm- For r > kn/2Nm, this probability is < /S^^/'^^^^ , 

Proof. Each site has k configs anchored at it, and, thus, the number 
of sites represented in the last r must be at least r^ = [r/k]. Now, select a 
subset of the n sites of size at least n/Nm consisting of sites none of which 
are in the m-neighborhood of each other. Assume < n/2Nm- The number 
of sites from the subset in the last is a hypergeometric random variable 
which, however, stochastically dominates a binomial random variable. Even 
assuming that all of the previous r^ — 1 sites selected have been from the 
subset size n/Nm, the probability the r^th is from the subset is > {n/Nm — 
r/k)/n > l/2Nm by assumption. Thus, the number of sites dominates a 
random variable X ~ Binom(rfc, l/2Nm). The conditional probability of any 
particular one of these sites being blocked is Pbiock < 1 — p + P{Am) < 1, 
using arguments from the last lemma. Thus, the probability of all of these 
being blocked is less than 

^{Pwock} = [1 - (1 - PUoc^)/2N^Y' < 

from which the first part of the result follows. Clearly, as r increases the 
probability can only decrease. □ 

In the next lemma we shall write E{r, N,m) instead of E{r,Qfyf ,m) since 
the bound does not depend on the positions of the sites in Qj^. 

Lemma 8. Assume that n is large enough that kn/2Nm = N/2Nm > 
n^/"^. Then the probability the last change in the RSA process on a set of 
n sites occurs at arrival N — r is < (3'^^^'' for r > n^/'^ . When r < n^/"^ the 
probability equals B'^~^^p{l — pY + Cr, where 



and e < B < e 



Proof. The first part follows directly from Lemma 7. Suppose now that 
r < n^/''. If the event E{r + 1, A^, m) occurs, then the last r + 1 arrivals are at 
sites not within m of each other. The conditional probability the last change 
is at kn — r is B'^~^^p{l —pY by Lemma 6. We then use Lemma 5 to get the 
bound on e^. □ 
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Proof of Theorem 1. The times between arrivals are independent of 
the order of arrival of the configs so that, if exactly j configs have arrived, 
the expected time until the next arrival is l/[kn — j). The expected time 
until the {kn — j)th arrival is thus X]f=j+i 1/^- It follows that the expected 
time until the last successful arrival is 



kn 1 kn—1 ( kn 



r=0 U=r+1 J r=[„i/7]_(_| U=r+1 J 

1/7 

where each ijr < f)'^ . As n — > oo, the last term tends to 0. Using Lemma 
8, the term in e breaks up into 

to U+i ^ J to'o^ 2^/iV t, 

and 

r=0 l[„i/7] M r=0 ViV 

Rearranging the last equation gives 

As n ^ oo, -B ^ 1, and this expression equals 

kn kn -. oo \j 

E7 + Mi-(i-P))-(i-P)t"^^^'-^ E 7+ E ^> 

«=1 j=[ni/7]+2 i=[nl/7]+2 

and the last terms tend to as n — > oo. □ 



3. The annihilating process. In RSA empty sites become filled as the 
process proceeds, whereas in annihilating processes, the sites begin occu- 
pied, and then some become unoccupied as time passes. In many cases a 
simple swapping of the roles of occupied and unoccupied turns one sort 
of process into the other. Annihilating processes have been considered in 
Sudbury [3] and Penrose and Sudbury [2]. Initially all sites of a graph are 
occupied. Between each pair of adjacent sites occur annihilation events at 
exponential rate 1. If one of the sites is unoccupied, nothing happens. If 
both are occupied, then one of the two sites, chosen at random, becomes 
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unoccupied. (By swapping the roles of "occupied" and "unoccupied," we 
have an RSA process in the sense of the previous sections, where a pair 00 
can be replaced by either 01 or 10.) If we start with n particles in a line and 
they evolve as an annihilating process, then at some time T„ the process will 
stop as all remaining particles are separated from each other. We solve for 
the distributions of the times in the sense that we find a generating function 
for them. Define 

(1) Fn{t)=P{Tn<t). 

Then it is clear that F(){t) = Fi{t) = 1. Further, if there are two separated 
blocks of m and n particles, then 

(t)Fn{t) = P(Tm, Tn < t). 

In a block of n particles, the end particles are annihilated at rate 1/2, and 
the interior particles are annihilated at rate 1, so the time until the first 
break is exponential with rate n — 1. Thus, 

(2) Fn{t)= y^Fr{t-s)Fr,.^_r{t-s)e-^^-^>ds, n>2. 

■^0 r=l 

Put Unit) = Fn{t)e'^\ then 

( n-l 

(3) n„(t)=e* / 'S^ Ur{t — s)un-i-r{t — s) ds , n>2. 
With U{x,t) = En=iUn{t)x'', we have 

oo n—1 

UrUn-r-lX^ = UiX + (uiMl + U2)x'^ + ••• = [/ + U'^ . 

n=2 r=l 

Thus, multiplying equation (3) by x" and adding from 2 to cxo, we obtain 

U-xe^ = xe^ f U{t - s) + U^{t - s) ds, 
Jo 

or 

t/e"* -x = x f U{s) + U^{s) ds, 
Jo 

and differentiating gives 

e-^^-Ue-' = x{U + U^). 
dt 



Putting V = l/U, we obtain 



— + (xe' + l)v = -xe\ 
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and using an integrating factor gives 

?; = -l + -e-* + c(rE)e-(^'^'+*\ 

X 

where c{x) is a function to be determined from the boundary conditions. 
Fi(0) = 1 and Fi(0) = 0, i > 1. Thus, C/(x,0) = x and 

(4) C/ = a;e*[l-x{e*-e-^'(^'-i)}]"\ 

We note that as t ^ cx), U — > xe^[l — xe*]"^ as would be expected since each 
Fi{t) 1. U{x,t) = J2T Fn{t){xe^)"-, so putting y = xe*, we obtain 



(5) 



Fiy,t) = J2Fnit)y'' = y[l-y + ye 



-t+y(e-*-l)]-l 



Equation (5) imphes 

Vfl - F (f))y- - y (e-V(l-j/))e(--*-^)^ 
Y "l-2/l + (e-V(l-y))e(-'-i)^" 

In order to find an asymptotic form for the mean time to stopping, we 
need to integrate this expression with respect to t. The term on the right- 
hand side is not directly integrable, but if we multiply by 1 + ye~* it is. We 
have 



Y^j {l + ye-%l-Fn{t))y^dt = - 



y 



In 1 + 



■ln{l-y} 



1-y 

(y + y2 + y3 + 



y2 y3 



SO that 



roo poo 2 1 1 

/ {l-Fn{t))dt+ e-*(l-F„_i(t))(it = l + - + - + ••• + -. 

JO Jo 1 6 n — I 

Equation (8) shows that Fn~i{t) — > as n — > cx), so that the second integral 
above — > 1. 

Theorem 9. The generating function for the duration of the annihila- 
tion process starting from n particles in a block is 



F{y,t) = J2Fn{t)y'' = y[l-y + ye 
1 

The mean time satisfies 



-i+y(e-'-l)]-l 



1 



1 1 

2 + 3 + -+n-l 
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